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Abstract 


An examination and brief review is made of the effects of quark-hadron 
transition induced fluctuations on big bang nucleosynthesis. It is shown 
that cosmologically critical densities in baryons are difficult to reconcile 
with observation, but the traditional baryon density constraints from 
homogeneous calculations might be loosened by as much as 50%, to 0.3 of 
critical density, and the limit on the number of neutrino flavors remains 
about < 4. To achieve baryon densities > 0.3 of critical density would 
require initial density contrasts R » 10 , whereas the simplest models for 

n 

the transition seem to restrict R to < 10 , 
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The possibility that effects due to the confinement of quarks in the 

1 2 

early Universe could create significant changes ' to the standard 

3 4 

homogeneous big bang nucleosynthesis results * has received a great deal of 
recent attention. In the standard homogeneous- Isotropic big bang 
nucleosynthesis calculation, a weak and nuclear reaction network is 
numerically followed for a uniform fluid, cosmologically expanding and 
cooling in the early universe to predict light element abundances. The 
success of such calculations is one of the central ingredients to the 
current overwhelming support found for the big bang model Itself. While the 
basic weak and nuclear reactions are measured in the lab to reasonable 
accuracy and are thus not seriously questioned, the assumption of a 
homogeneous- isotopic fluid has been questioned many times (cf. ref. 3). 
Recent work on the quark-hadron transition has given a physically derived 
motivation to such questioning. In particular the transition from the early 
"quark- soup” to normal hadronic nuclear matter should take place at T > 100 
MeV at just prior to the nucleosynthesis epoch, at T < 1 MeV. Witten and 
others^ had noted that if the quark-hadron transition is a first order 
phase transition then density fluctuations would naturally result. The 
possible effects of these fluctuations on big bang nucleosynthesis 
calculations is the reason for the current excitement. 

The purpose of this paper is to briefly review the previous quark- 
hadron inspired results and compare them with the traditional homogeneous 
results and then to present a new set of calculations which explicitly show 
the sensitivity of the resultant light element abundances to the parameters 
of the quark-hadron transition. We will show that even If the transition is 
first order, the result is unlikely to signif icnatly alter the key 


-4- 


predictions from homogeneous nucleosynthesis as long as one continues to 

require agreement with the observed light element abundances, particularly 
7 4 

Li and He, The persistence of the nucleosynthesis conclusions despite the 

addition of new initial conditions with several additional parameters shows 

the robustness of big bang nucleosynthesis. 

Traditional big bang nucleosynthesis had become one of the cornerstones 

of big bang cosmology because of its remarkable agreement with light element 

abundance observations, spanning a dynamical range of over 9 orders of 

magnitude in its predictive powers. This success, coupled with its 

6 7 3 

prediction of the number of neutrino families ’ ’ is an important 

vindication of the "particle physics connection" in the study of the early 

universe. Furthermore, standard big bang nucleosynthesis arguments using 

deuterium 8 ^ and later helium-3 3 ^ and lithium 3,9 ^ constrain the density, 

of normal matter, baryons^ , in units of the critical density to ^ - 0.1. 

More precisely, the ratio of baryons to photons, n^/n - tj is 
3 9) 

constrained ’ to 


3 x 10' 10 < n < 4 x 10' 10 


( 1 ) 


for current population II stellar lithium abundances and current limits on D 
3 * 

and He , The fact that ~ 1 is excluded is one of the prime driving 
forces behind the current searches for non-baryonic dark matter^^ . 

* -10 7 

ThjgUpper limit of rj < 6 x 10 in ref. 3 from Li was reduced to i) <’5 x 

10 in ref. 9 using newer lithium rates. The latest rates^" ' yield the 

upper limit In eq. (1). 
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1 2) 

With so much at stake , the initial claims * that a quark-hadron 

transition inspired model could yield an - 1 Universe compatible with 

light element abundances created tremendous interest. Some preliminary 

lattice gauge calculations implied that the quark-hadron transition may 

indeed be a first order phase transition. Applegate et al^ noted that due 

to the proton's electric charge there is preferential diffusion of neutrons 

versus protons out of the high density fluctuations produced by such a 

quark-hadron transition. This could lead to big bang nucleosynthesis 

occurring under conditions with both inhomogeneities and variable 

neutron/proton, n/p, ratios. The result is that the nucleosynthesis in the 

high density regions occurs with a low n/p ratio while the low density 

region has a high n/p. Regions with n/p > 1 have qualitatively different 

nucleosynthesis then standard homogeneous nucleosynthesis (where n/p - 1/7). 

If n/p > 1, the number of protons rather then neutrons becomes the 

4 

constraining parameter on the reaction network flow towards He. 

1 2 ) 

In the first round of calculations * 7 these groups claimed that such 

mixed conditions might allow Ci^ ~ 1 while fitting the observed primordial 

4 3 7 7 

abundances of He, D, He but with an overproduction of Li. Since Li Is 

the most recent of the cosmological abundance constraints and has a 

different observed abundance In population I stars versus the traditionally 

12 ) 

more primitive population II stars ' some argued that perhaps some special 
depletion process might have occurred to reduce the excess ^Li. Reeves and 
Audouze et al 7 each argued against such processes and tried to turn the 
argument around and use the lithium abundances to constrain properties of 
the quark-hadron transition. 
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On this basis , Reeves concluded that the ^Li abundances required that 

the ratio R, of baryon densities in the high to low density regions satisfy 

R <2-4. These limits in principle imply constraints on the transition 

temperature T £ 150 MeV. The limit on T is however based on naive 
c c 

assumptions made in estimating the density contrast as a function of the 

5 2 ) 

transition temperature T 1 . . The main ingredient neglected was the 
interactions in the hadron phase (indeed without these one would conclude 
the existence of a high temperature hadron phase) . When the effects of the 
finite size of hadrons due to repulsive interactions are included 1 ^ one 
finds that for a first order transition R > 7 for all values of T This 

- c * 

means that possible constraints from nucleosynthesis must be on the more 
detailed aspects of the phase transition. One should also note that the 
baryon density contrast across the phase boundary during the transition does 
not necessarily translate directly into the density contrast remaining after 
the transition^^ . 

At first it appeared that if the lithium constraint could be surmounted 
then the constraints of standard big bang nucleosynthesis might 
disintegrate. Although the number of parameters needed to fit the light 
elements was somewhat larger for the non-standard models, nonetheless a non- 
trivial loophole appeared to be forming. To further stimulate the flow 
through the loophole, Malaney and Fowler 16 ^ showed that in addition to 
looking at the diffusion of neutrons out of high density regions one must 
also look at the subsequent effect of neutrons diffusing back into the high 
density regions as free neutrons are depleted at a much slower rate in the 
low density regions in nucleosynthesis . (The initial calculations treated 
the two regions separately.) Malaney and Fowler argued that for certain 
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phase transition parameter values, (eg. nucleation site separations - 10m at 

the time of the transition) this back diffusion could destroy much of the 

7 7 7 4 

excess lithium produced as Be via Be(np) Li(pa) He in the high density 

17 18 19 ) 

regions. However, it has been recently argued * * } that in detailed 

diffusion models, the back diffusion not only affects ^Li but also the other 

4 

light nuclei as well. Those calculations found that for Q^~l, He is also 
overproduced (although it does go to a minimum for similar parameter values 
as does the lithium) . 

4 7 

One can understand why these models tend to overproduce He and Li by 

remembering that in standard homogeneous big bang nucleosynthesis, high 

baryon densities lead to excesses in these nuclei. As back diffusion evens 

out the effects of the initial fluctuation the averaged result should 

approach the homogeneous value. Furthermore, any narrow range of 

parameters, such as those which yield relatively low lithium and helium, are 

unrealistic since in any realistic phase transition there is a distribution 

of parameter values (distribution of nucleation sites, separations, density 

20 ) 

fluctuations etc.). Therefore narrow minima are washed out which would 
7 4 

bring the Li and He values back up to excessive levels for parameter 
values with 0^-1, We stress this point since diffusive effects are only 
important in lowering the nuclear abundance in a narrow window of parameter 
space . 

After the above review of the current situation and the new apparent 
difficulties in making ~ 1, we have decided to address the quark-hadron 
transition with a more traditional approach. Namely, instead of setting 
- 1 and seeing what excesses may or may not occur, let us believe the light 
element abundance observations and see how the traditional big bang 
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nucleosynthesis constraints might vary as quark-hadron transition parameters 
are explored, (This is similar to the approaches of Reeves and Audouze et 

13 ) 

al however we are using the more detailed dynamical code of Kurki-Suonio 

lg\ 

and Matzner which explicitely includes multizone forward and backward 

diffusion). Indeed, one might worry that because R £ 7 for all values of 

T c . the allowed set of parameters In standard big bang nucleosynthesis might 

be altered (e.g. the range In iy) . We will therefore test the standard model 

parameters in the presence of baryon inhomogeneities. 

In these calculations we did not explore the exciting possibility^^ 

that quark-hadron fluctuations might enable big bang nucleosynthesis to make 

elements heavier than ^Li which are blocked in the conventional model. If 

such synthesis is possible for the allowed parameter space that fits the 

light element abundances this would be very exciting and might explain some 

abundance patterns In metal-poor stars and provide an Independent test of 

whether or not the transition was indeed first order. 

We have also not explored the remaining fundamental physics questions 

about the transition itself. Is it a first order phase transition? What is 

the relationship between nucleation sites, density fluctuations, etc., and 

the fundamental QCD parameter We have also followed the previous 

calculations and assumed basically Isothermal fluctuations, however 

differential temperature diffusion should be explored. 

We follow the parameterization of Kurki-Suonio and Matzner which 

treats the transition in a very phenomenological manner. (For the 

relationship of these parameters to certain bag models see Alcock 

2 ) 

et al .) Thus our aim in this paper is not to make specific statements 
about the physics of the quark-hadron transition (although some inferences 
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might be made) but instead to see what effects the transition might have on 
the traditional big bang nucleosynthesis constrained quantities, in 
particular on f|, To this end, we will use the results of Kurki-Suonio and 
Matzner^^ for the calculated abundances of D, ^He, Sie and ^Li in a 
nucleosynthesis model in the presence of baryon inhomogeneities with 
diffusion taking place before and during nucleosynthesis. Because the 
details of the quark-hadron transition are largely unknown, we explore a 
parameter space to find the largest possible set of primordial abundances. 

The phenomenological parameters we explored which can affect 
nucleosynthesis are the following: 

(1) The average baryon to photon ratio, »?; 

(2) The average density contrast R 

(3) The average distance, scale of the inhomogeneities, i and 

(4) The average volume fraction of the high density regions, f . 

(Note that only r? is a parameter in the homogeneous case.) Furthermore the 

geometry of the high density regions can also have an effect. We consider 

planar, spherical and cylindrical geometries. In this paper we did not 

consider fractal -like boundaries which might also result in such transitions 

and could further enhance surface diffusion effects. Our results are 

displayed for a baryon density contrast between the high and low density 

phases, R - 100. Increasing (decreasing) R, brings the resultant abundances 

18 \ 

further from (closer to) the homogeneous results , For example, reducing 
the contrast to R - 10 depending on the volume fraction involved, reduces 
the deviation from the homogeneous results to about one -half and to about 
one-quarter for R - 6 (when f - ^) . We also consider a range - 1/4 - 
1/64 for the volume fraction of the high density region. Specifically, our 
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results use data from the following choices of parameters - 1 / 4 , f - 1/8 
and f^ - 1/16 for planar geometries and f^ - 1/8 and f^ - 1/64 for spherical 
geometries. These were chosen so as to minimize and maximize the elemental 
abundances. As f goes to 1 or 0, the results approach the homogeneous 
results, (For larger R, smaller f would have to be considered.) The 
distance scale £ is given in meters at 100 MeV after the phase transition 
between the centers of high and low density regions. 

For a given value of t) and £, we have varied f and the geometry so as to 
find a maximal range for the calculated abundances. We will find that only 
for a limited range in £ and »? are the derived abundances in agreement with 
observational determinations. 

In the figure , we show the allowed region in the £ - jj plane from the 
constraints given by the abundances of D, 3 He, ^He and 7 Li. The 

observational constraints we use are the following 4 ^: D/H 10 ‘ 5 * by number, 

3-4 / 

(D + He)/H < 10 by number, 0.224 i Y 4Hg < 0.254, where Y„ He is the *He 

abundance by mass and 7 Li/H < 2 x 10 ^ by number for population II and 

7 .9 

Li/H < 2 x 10 by number for population I. For standard big bang 

nucleosynthesis (i - 0 ) the bounds on r} may be read from the bottom of the 

figure; they are the results giving rise to eq. ( 1 ). 

The calculated abundances are for a neutron half-life of r - 10 35 

n 

min. The weak n <--> p rates are obtained by numerical integration, and 
multiplied with a Coulomb correction factor 22 ^ 0.98. Additional small 

23} 

corrections calculated by Dicus et al. ' are represented by subtracting 

4 

0.001 from all He mass fractions. The strong reaction rates used are from 

the recent compilation by Caughlan and Fowler^^ . The new rates for 

2 32 33 44 7 7 7 

H(d,n) He, H(d,p) H, He(d,p) He, He(t, 7 ) Li, and Be(n,p) Li lead to a 
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higher estimate for produced ^Li further narrowing the range of *7 allowed by 

7 7 4 

population II Li. Rates for (n,7) reactions and Be(n,a) He not included 

24) 

in this compilation are those used by Schramm and Wagoner , except the 

7 8 25 ) 

newer estimate for Li(n,7> Li by Malaney and Fowler } is used. Abundances 
for A > 7 isotopes are not calculated but their maximum effect on A < 7 
isotopes was controlled by including the reactions leading to A > 7 as 
sinks. For the density range discussed here the effect of these sinks on 

final ^Li was at most a few per cent. (Except that in the f - 1/64, ij - 7 

-9 7 

x 10 case, where the high density region had the highest density, Li came 

15-45% lower with sinks than without sinks. A full network would give a 

result in between. Since these yields were an order of magnitude above the 

population I upper limit, this inaccuracy does not affect the results 

reported here.) Because reactions occur in thin layers near the original 

high/low density boundary 4 * * * * 9 , a fairly fine zoning was necessary for 

accurate results (in most cases 64 zones was found to be sufficient, 

20 ) 

compared with only 8 zones for the Livermore group and 2 for the Tokyo 
19), 

group ) . 

When we compare these to the observational constraints, we find the 

contours shown in the figure. Consider for example the contour found from 

4 

< 0.254. As the inhomogeneity is turned on, the He abundance 

4 

increases. The rise in He for small i has a straightforward explanation. 
Because of the small distance all the neutrons can find their way to a 

4 

high density to react prior to their decay. He is raised in the high 

density regions. Since the computations are for fixed averaged baryon 

density, the result is that nucleosynthesis occurs in overdense regions 
4 

giving enhanced He. For very small i, the protons as well as the neutrons 


4 
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diffuse, giving the uniform baryon density result when nucleosynthesis 

begins. For an optimal value of i in the range 10-100, one sees the 

or l-8i riS l effect claimed in refs. 1 and 2. However, as claimed there, the 

back diffusion does not allow the drop in Y 4 ^ e to be as pronounced and only 

a modest increase in the limit on tj based on Y 4 ^ is seen. For larger 

values of f, diffusion becomes irrelevant and one has strictly an 

inhomogeneous nucleosynthesis model and one finds a larger Ste 
3 27 6 ) 

abundance ’ 1 (and hence a tighter constraint on 17). 

3 

For the cases of D and D + He, aside from a slight decrease in D (for 
relatively low r;) both D and D + 3 He increase with i. This shifts allowed 
the values of jj to a higher range. For i — 0, standard nucleosynthesis, D 
and D + 3 He require 3 x 10‘ 10 £ rj < 10 x 10" 10 whereas, for i 2: 100 this 
range moves up to 4 . 5 x 10' 10 < ij < 30 x 10‘ 10 . The dip in D for i - 10 
allows a drop in the bound in r) , r/ > 2.2 x 10" 10 . 

The 7 Li abundances, as has been known all along in this type of 
investigation, rise with increasing £\ the effect of which is to decrease 
the allowed range for t] . In the case of the population II 7 Li abundances, 
we see rather dramatically the constraint £ < 150, for any value of rj . When 
t ~ 10, we find (using the D + He abundances for the lower limit) 2.2 x 
10 i rj < 3 x 10 There is a gap which excludes values of i from 30- 

100. For i between 100 and 150, we have 4 x 10' 10 £ *7 s 7 x 10* 10 , limits 
which are comparable or tighter than the standard nucleosynthesis bounds. 

The bound from 7 Li JI , t) < 7 x 10' 10 for £ ~ 100 is evidence of the Malaney- 
Fowler^^ effect, a maximization of back diffusion destruction of 7 Li 
(actually Be, which produces Li by e -capture). The 7 Li abundance is 
determined late in nucleosynthesis, when neutron abundance is very low. For 
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- y 

17 < 10 the neutron fraction in the low- density region is diminishing 

slowly enough that sufficient neutrons are available to diffuse into the 

high-density regions and destroy most of the 7 Be there. This effect is very 

sensitive to the distance scale. If distances are too short, diffusion 

depletes the neutron resevoir too early. If distances are too long, 

diffusion will not be efficient. For higher densities, e.g. rj — 7 x 10"^ 

with R - 100, even the low-density region is too dense for a sufficient 

number of neutrons to survive long enough to have a dramatic effect. For 

the population I abundances , though we do not find a limit on i , the bounds 

on ri are again comparable to the standard results. In either case, r) < 20 x 
-10 

10 or (^ < 0,3 remain upper bounds for all values of the parameters 
considered. We conclude once more that the Universe can not be closed by 
baryons. (The lower bound of drops by only - 25%; thus still being 
greater than 0 in visible matter.) 

Diffusion effects on nucleosynthesis could be stronger if the density 

contrast were much higher than R — 100, Because the details of the 

confinement transition are poorly understood, it is difficult to make a 

convincing calculation of R from first principles. One approach 5 ^ has been 

widely used; namely calculate R assuming chemical equilibrium during the 

phase transition. With this assumption (and only with this assumption) can 

one calculate unambiguously the density contrast. In this case, it was 
14) 

shown 7 that for T c > 100 MeV, R < 100. 

Kurki-Suonio ^ considered possibilities for the evolution of baryon 
number fluctuations assuming that the equilibrium ratio is maintained at the 
phase boundary but only extends a diffusion length from the boundary. 
Depending on the distance scales of nucleation, coalescence, and diffusion, 
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he obtained various possibilities, with the most probable being to form 

final density contrast R - (w /w )R , where (w u /w_) is the ratio of 

rl eq n 

enthalpy densities of the two phases, which is less than 1, and R is the 

eq 

equilibrium baryon density ratio. To obtain significant inhomogene I ties 
with R much larger that R would seem to require extremely efficient baryon 
transport in the quark phase, the more likely outcome being that the final 
inhomogeneity involves only an Insignificant fraction of the total baryon 
number. 

In the models of Ref. 18, the dependence on R Is rather weak and 
results for R - 1000 do not appear very different. In. Ref. 19, R - 10^ - 

4 

10 , was claimed to allow 0^ - 1 if h Q < 0.5. However this possibility is 

achieved only for an extremely narrow range in the parameter f^. (We remind 

the reader of our previous comment with regard to results which are valid 

only in narrow windows.) This conclusion is based on a two-zone calculation 

(in contrast to the 64-zone calculation in ref. 18) and uses constraints 

Y 4He <0.26 and 7 Li/H < 10 9 (we assume that Y 4Hg < 0.254 and 7 Li/H < 2 x 

10 10 for population I and 7 Li/H < 2 x 10* 9 for population II). The 

homogeneous value of Y 4Hg in ref 19 also falls short by about 0.005 of the 

homogeneous calculations used here and in ref. 3. On this basis, we do not 

feel that there is any real disagreement between those results and the ones 

28 ) 

quoted here. Mathews et. al . have studied the effect of extreme density 

5 2 3 

contrasts R - 10 and report that with suitable parameter values H, He, 

4 

He can be brought to simultaneous agreement with observations. From 

comparisons with the work of other groups, it would seem that deviations 

3 

from our conclusions only begin to occur for R » 10 which we consider 


unrealistic . 
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It Is also interesting to note that in addition to altering the bounds 

on t} , baryon inhomogeneities and neutron diffusion could In general alter 

the limits on the number of neutrino flavors from nucleosynthesis. The 

current limit of N <4.2 based on Y 4l , < 0,254, r >10.2 min and n > 3 x 

v 4 He n 

10 Depending on the value of rj and i, the limit could increase or 

decrease. For example, at rj - 4 x 10" 10 and £ - 100, < 4,6 and i ~ 10 

for the same value of rj } N <3.9, while for q - 2.2 x 10 ^ with i - 10, N 

v v 

<4.3. Thus again, we find only minor fluctuations from the traditional 
conclusion. 

Although from nucleosynthesis abundances we can not calculate a limit 

to R or the possibly related parameter T , the limit on £ is an Interesting 

constraint (though not a terribly strong one) on the quark hadron 

transition. The distance scale £ has been estimated in terms of transition 

parameters such as the transition temperature T , the surface tension 

associated with the fluctuations, a, and the latent heat of the transition, 

L. Assuming L - 15 T Fuller et al 29) find i - (4 x 10 4 ) (a/MeV 3 ) 3/2 
-13/2 

(T^/MeV) . (We have here corrected for the error in the numerical 

factor in the approximate solution for the supercooling parameters in Fuller 

et. al, which was too large by a factor - 4, making their distance scale 

30) 

estimates 50 times too large , The surface tension a has been 
31) 1/3 

estimated o' <70 MeV so that for > 100 MeV we expect that i < 1, 
(and note the strong temperature dependence) well below our nucleosynthesis 
bound of £ < 150. 

In conclusion, we find that for reasonable values of the baryon density 
contrast R < 100, it remains possible to be consistent with observational 
determinations of the light element abundance (including population II ^Li) 
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if the mean separation of the fluctuations is i < 150. In addition the 
standard nucleosynthesis constraints on r} and remain largely intact. For 
i “ 10, may be as low as 2.2 x 10* 10 (but less than 3 x 10" 10 ). The upper 
bound on r\ is r\ < 7 x 10’ 10 for all values of i for population II 7 Li 
abundances. This upper limit is increased to r\ < 20 x 10’^° for the 
population I 7 Li abundances. In all cases we find - 1 still excluded by 
big bang nucleosynthesis . 
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Figure Caption 

Allowed regions in the i-r? plane from the observational constraints on 
3 4 7 

D, He, He and Li (from both population I and II stars). The area 
outlined by bold lines are the only regions consistent with all 


observations , 
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